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! We propose a new thermalization scenario for heavy ion collisions which at suffi- 

' ciently high energies implies the phase transition to the quark-gluon plasma. The 

, key ingredient of our approach is the Hawking-Unruh effect: an observer moving 

, with an acceleration a experiences the influence of a thermal bath with an effective 

cn ■ 

CS| ' temperature T = a/27r, similar to the one present in the vicinity of a black hole 

■ horizon. For electric charges moving in external electromagnetic fields of realistic 

■ 

i strength, the resulting temperature appears too small to be detected. However for 

Q^' partons in strong color fields the effect should be observable: in the Color Glass 

Condensate picture, the strength of the color-electric field is E ~ Q'i/g {Qs is the 

• • . saturation scale, and g is the strong coupling), the typical acceleration is a ^ Qs, 

> ■ 

^ I and the heat bath temperature is T = Qs/2tt ~ 200 MeV. In nuclear collisions at 

' sufficiently high energies the effect can induce a rapid thermalization over the time 

period of r ~ 27r/Qs ~ 1 fm accompanied by phase transitions. We consider a spe- 
cific example of chiral symmetry restoration induced by a rapid deceleration of the 
colliding nuclei. We argue that parton saturation in the initial nuclear wave functions 
is a necessary pre-condition for the formation of quark-gluon plasma. We discuss 
the implications of our "black hole thermalization" scenario for various observables 
in relativistic heavy ion collisions. 
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I. INTRODUCTION 

In 1974 Hawking Q demonstrated that black holes evaporate by quantum pair produc- 
tion, and behave as if they have an effective temperature of 

= h 

/TT 

where k = {AGM)~^ is the acceleration of gravity at the surface of a black hole of mass 
M; G is the Newton constant. The thermal character of the black hole radiation stems 
from the presence of the event horizon, which hides the interior of the black hole from an 
outside observer. The rate of pair production in the gravitational background of a black 
hole can be evaluated by considering the tunneling through the event horizon. Parikh and 
Wilczek showed that the imaginary part of the action for this classically forbidden process 
corresponds to the exponent of the Boltzmann factor describing the thermal emission*. 

Unruh ^ has found that a similar effect arises in a uniformly accelerated frame, where 
an observer detects an apparent thermal radiation with the temperature 

Tu = ^; (2) 

(a is the acceleration). The event horizon in this case emerges due to the existence of 
causally disconnected regions of space-time jj] , conveniently described by using the Rindler 
coordinates. 

The effects associated with a heat bath of temperature usually are not easy to de- 
tect because of the smallness of the acceleration a in realistic experimental conditions. For 
example, for the acceleration of gravity on the surface of Earth g ~ 9.8 m s^^ the corre- 
sponding temperature is only T ~ 4 x 10~^° K. A much larger accelerations can be achieved 
in electromagnetic fields, and Bell and Leinaas considered the possible manifestations 
of the Hawking-Unruh effect in particle accelerators. They argued that the presence of 
an apparent heat bath can cause beam depolarization. Indeed, if the energies of spin-up 

and spin-down states of a particle in the magnetic field of an accelerator differ by 
AE = Ei^ —Et^, the Hawking-Unruh effect would lead to the thermal ratio of the occupation 



Conservation laws also imply a non-thermal correction to the emission rate ^ , possibly causing a leakage 
of information from the black hole. 
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probabilities 




AE 



) 



(3) 



where the Unruh temperature Q is determined by the particle acceleration. According to 
(121), a pure polarization state of a particle in the accelerator is inevitably diluted by the 
acceleration. 

If the energy spectrum of an accelerated observer is continuous, as is the case for a particle 
of mass m with a transverse (with respect to the direction of acceleration) momentum pfi 
a straightforward extension of leads to a thermal distribution in the "transverse mass" 
niT = \Jni? + p^: 



An important example is provided by the dynamics of charged particles in external electric 
fields. Since a particle of transverse mass rriT and charge e in an external electric field of 
strength E moves with an acceleration a = eE/mT, the thermal distribution ^ with the 
temperature Tu = a/2n is given by 



which differs from the classical Schwinger result |6[ for the momentum distribution of par- 
ticles produced from the vacuum only by a factor of two in the exponent. We will establish 
in Section inithat this interpretation of Schwinger's result indeed holds, and will show how a 
more careful treatment avoids the factor of 2 discrepancy. Note that the result (0) cannot be 
expanded in powers of the field strength E, and thus cannot be reproduced in perturbation 
theory. This fact has a geometrical interpretation in the accelerated frame: the Bogoliubov 
transformation relating the particle creation and annihilation operators in Minkowski and 
Rindler spaces describe a rearrangement of the vacuum structure which cannot be captured 
by perturbative series. 

Perhaps the largest accelerations accessible to experiment at present are achieved in the 
collisions of relativistic heavy ions. Indeed, at RHIC accelerator at BNL, the heavy ions 
collide with c.m.s. momenta of 100 GeV per nucleon, and the strong interactions of their 
parton constituents lead to the production of a sizable fraction of final state partons at 
rest in the c.m.s. frame. This happens over a very short time At < 1 fm, so a typical 
deceleration is a ~ (At)^^ 200 MeV (we will find later that due to the presence of strong 




(4) 
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color fields in the initial state described by parton saturation the achieved deceleration is 
even higher, a ~ 1 GeV). Such deceleration should lead to observable effects, which include 
an apparently fast thermalization of the produced partonic state. Indeed, the experimental 
data from RHIC 

mm 

and their analysis in terms of the hydrodynamical approach 
[ll . 12, 12, 14, 3| suggest that thermalization occurs over a time of about 1 fm, much faster 
than is expected on the basis of perturbative calculations Q|. 

In this paper we will argue that the Hawking-Unruh effect indeed can cause rapid ther- 
malization in relativistic heavy ion collisions. If we describe the initi al p arton wave function 
of the nucleus in the Color Glass Condensate picture 

HHQ, HHQ, HQS 

where the strength of the color-electric field is E ~ Q^/g {Qs is the saturation scale, and 
g is the strong coupling), the typical acceleration achieved in a collision is a ~ Qs, and the 
resulting heat bath temperature is T = Qsf^Tr ~ 200 MeV. As will be discussed below, in 
nuclear collisions at sufficiently high energies and/or sufficiently large atomic number, the 
effect can induce a rapid thermalization over the time period of r ~ 27r/(5s accompanied by 
phase transitions. 

The paper is organized as follows: In section |n] we consider the case of particles moving in 
strong external fields of various kind, and discuss the correspondence between the Schwinger 
and Hawking-Unruh effects. We show in particular that a pulse of a strong (color)electric 
field produces thermally distributed charged particles. In section IIIII we formulate an ap- 
proach to relativistic heavy ion collisions which takes into account the Hawking-Unruh effect. 
Soft partons in this picture are produced by tunneling through the event horizon in Rindler 
space, and emerge thermally distributed with temperature T = Qs/^tc ~ 200 MeV. We de- 
rive the notion of a critical acceleration Oc within the dual string picture, and show that it is 
related to the Hagedorn temperature Tuagedom by the relation Tnagedom = ac/27r. We argue 
that parton saturation in the initial wave function is a necessary condition of deconfinement 
in the final state, since it allows to create decelerations larger than Oc, and thus to produce 
a thermal medium with a temperature T > THagedom- As an explicit example of a phase 
transition caused by a strong deceleration a ~ Q^, we discuss in section HVBI the restoration 
of chiral symmetry in the framework of Nambu-Jona-Lasinio model formulated in Rindler 
space. We then discuss the implications of our scenario for various observables in relativis- 
tic heavy ion collisions. We conclude the paper with the Summary, where we discuss the 
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limitations of our approach, its relation to other developments, and suggest some directions 
for future studies. 



II. PAIR PRODUCTION IN STRONG EXTERNAL FIELDS: TUNNELING 

THROUGH THE EVENT HORIZON 

A. Schwinger pair production as a Hawking Unruh phenomenon 

In a classic 1951 paper Q], Schwinger considered the dynamics of QED in the strong field 
domain. He constructed an effective action describing the coupling of electromagnetic fields 
to charged particles; in an electric field the action acquired an imaginary part corresponding 
to the pair creation. For future discussion it is instructive to rederive here the result of 
Schwinger in the Euclidean time formalism. Let us consider the action of a charged particle 
of mass m in a constant external electric field E: 

S = J { — mds — e(fdt) , (6) 

where cp is the electric potential. In a constant electric field E the electric potential is 

if = —Ex modulo an additive constant. Equations of motion of this particle are 

dpx dp± , . 

— — = eE, —— = 0. (7) 
dt dt ^ ' 

The velocity v of the particle is v = dx/dt = p/ £\a-ai where the kinetic energy of the particle 

is = rn? + fP' . With the initial conditions p(0) = 0, x(0) = equations of motion ((7j) 

can be integrated yielding the trajectory 

vit) = , (8a) 



x{t) = a-^ ( Vl + aH^ - 1) , (8b) 

where we can identify the constant a with acceleration of the particle: 

d V eE 
a = — , = — ; (9) 

dt Vl - m ' ^ ^ 

the trajectory ()8b|) is shown by the upper curve in Fig.l. 

At t < the particle is coming from x oo gradually losing its kinetic energy in favor 

of the potential one and finally stopping at x = at t = 0. This is the turning point of the 
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classical trajectory. At t > the particle accelerates again to the velocity v ^ 1 at x oo. 
The action associated with the particle moving along the trajectory ()8b|) is given by ®. 
Using ds"^ = (1 — dt"^ and substituting the trajectory (|Hajl and (j8b|l into ^ we get 

S{t) = J dt{-m^l - v{tY + eEx{t)) 

= arcsinh(ar) + — - ( a r (Vl + a? t"^ — 2) + arcsinh(ar)) + const. (10) 

a 2 a'' \ / 



In classical mechanics the equations (jHaj) . (j8bjl completely specify the motion of a uniformly 
accelerating particle moving under the influence of a constant force F = — eVy?- In contrast, 
in quantum theory the particle has a finite probability to be found under the potential 
barrier V{x) = eEx in the classically forbidden region. Mathematically, it comes about 
since the action ()10|) being an analytic function of r has an imaginary part 



lmS{r) = ^-'-^ = ^. (11) 
^ ' a 2a^ 2eE ^ ^ 



The imaginary part of the action ()10j) corresponds to the motion of a particle in Euclidean 
time. Substituting t —He in (|8bj) we find the Euclidean trajectory 



x{tE) = a-'i^l - aHl - 1). (12) 

Note that unlike in Minkowski space the Euclidean trajectory is bouncing between the two 
identical points Xa = —a~^ at tE^a = — and x^ = — a^^ at tE^b = and the turning 
point = at t^; a = - see the lower curve in Fig.l. Using (fTUI) we can find the Euclidean 
action between the points a and b; it is given by SE{x{tE)) = T^im? /2eE. It is well known 
that a quasi-classical exponent describing the decay of a metastable state is given by the 
Euclidean action of the bouncing solution (see e. g. {2^), (fT^ . The rate of tunneling under 
the potential barrier in the quasi-classical approximation is thus given by 

Vy^m - e~''-' = e-^. (13) 

Equation (fT!?|) gives the probability to produce a particle and its antiparticle (each of mass 
m) out of the vacuum by a constant electric field E. The ratio of the probabilities to produce 
states of masses M and m is then 

Ty^M _ ^(A/2-^2) 

= e . (14) 

J- V^m 
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ows a dual interpretation in terms of both Unruh and Schwinger 



The relation (HJl al 

effects (see e.g. js?! 2^ 29| and references therein). First, consider a detector with quantum 
levels m and M moving with a constant acceleration in the constant electric field. Each level 
is accelerated differently, however if the splitting is not large, M — m <^ m we can introduce 



the average acceleration of the detector 

2eE 



(15) 



M + m 

Substituting (fT3j) into (fT^ we arrive at 

= e ^ . (16) 

J- V^m 

This expression is reminiscent of the Boltzmann weight in a heat bath with an effective 
temperature 1^: T = a/2TT. It implies that the detector is effectively immersed in a photon 
heat bath at temperature T ^ eE/nm. This is the renown Unruh effect j^. 

Schwinger mechanism traditionally is viewed as a process of creation of virtual particle- 
antiparticle pairs from vacuum. In the case of electrodynamics this corresponds to the 
tunneling of an electron from the Dirac sea |^. In this process the electron energy 
changes from £_ to £+, where 

e± = ±A/p2(x) + m2 + \e\Ex. (17) 

When the quasiclassical action is evaluated along the true trajectory it can be viewed as a 
function of the final coordinate (or time) of the moving particle. In this case its imaginary 
part is given by ImS* = J dx \p{x)\ where the integral must be evaluated in the region of 
imaginary p between turning points a, h such that p{xafi) = 0, i. e. Xa^ = (e ± m)l\e\E. A 
simple calculation yields the result for the tunneling probability which coincides with (fT^ . 

Another interpretation of (fT^ can be given if we recall that the uniformly accelerated 
detector in Minkowski space is equivalent to the inertial detector in the Rindler space. The 
vacuum in Minkowski space is related to the vacuum in the Rindler space by a non-trivial 
Bogoliubov transformation which shows that the Rindler vacuum is populated with the 
thermal radiation of temperature T = 27r/a in accordance with In general, according 

to the Equivalence Principle of General Relativity, equations of motion in an accelerated 
frame are equivalent to a motion in a certain gravitational field. The thermal spectrum 
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FIG. 1: Trajectories of an accelerated particle in Minkowski and Euclidean space; the spatial 
coordinate x and the time t are plotted in units of the inverse acceleration, . The upper curve 
represents the Minkowski trajectory of a particle under the influence of a constant force F: the 
particle approaches from x = +oo, slows down and stops at a; = 0, i = 0, then turns around and 
accelerates again towards x = +oo. The trajectory is a hyperbola in Minkowski space, or a line 
with constant value of the coordinate = x"^ — = in Rindler space. The lower curve is 
the instanton-like Euclidean trajectory of the particle under the potential barrier V = —Fx: the 
particle bounces between x = — 1/a and x = with a period of motion (3 = In ja. The inverse of 
/3 has a meaning of an effective Unruh temperature Tu = 

appears then as a result of tunneling through the event horizon. The most famous example 
of such a tunneling process in quasi-classical quantum gravity is Hawking radiation Q where 
the thermal quanta are emitted from the black hole horizon. We would like to note that the 
methods of General Relativity have been applied recently to the studies of different aspects 



The method we used for derivation of the Schwinger formula yields another important 
well-known result: constant magnetic field cannot produce pairs from vacuum. Indeed, an 
electron in this case moves with a constant absolute value of the acceleration on a circular 
orbit. In such a motion = const and hence ImS* = 0. 



of QCD by several authors |32l. l33l l34 . 
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B. Pair production by a pulse of electric field 

In the situations we consider in this paper the fields can be assumed approximately spa- 
tially homogeneous. However, the interaction time is short, so the fields are time-dependent. 
Let us therefore use the method described above to calculate the rate of pair production by 
a spatially homogeneous but time-dependent electric field 

E = Excoshr^{t/i). (18) 

Our choice for the functional time dependence of the field pulse is motivated by the fact 
that Klein-Gordon and Dirac equations in this case are exactly integrable (see e.g. 
and we can test the accuracy of our quasi-classical result. 

The equations of motion of an electric charge e in this case are given by 

^^ = ^dkNy ^^ = °- ^''^ 

With the initial condition p(0) = they are integrated to give 

p^{t) = e^t tanh(t/t) , py{t) = 0. (20) 



The velocity of the particle is given hy v = p/S]^i^ = pj ^vn? + pi- Therefore, according to 
(fTn|l . the kinetic part of the action along the trajectory is 



^ , , -m^t , Vm^ + e2 tanh(t/t) , , 

Sk{t) = arctanh ^ ' L. (21) 



The imaginary part of the action (j2T|l is 



Im SKir) = , """"'^ . (22) 

The contribution of the dynamical part to the action is given by the second term in the 
right-hand-side of (jHI). To find an explicit form of the trajectory x{t) let us introduce a new 
variable ^ = sinh(t/t). Then, 

x(0 = eEe ^ "^^^ 



1 + ^2 + ^2 ^^2 + g2 ^2 ^2) 



+ e2 E^ t2 P V m2 + e2 ^2 t2 



Ini x/TTe + \ e + -n^^^-TT^ } + const. (23) 
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Correspondingly, the action is 
Sj{t) = eE J x{t)dt -- 



j dt\n\ cosh(t/t) + y sinh2(t/t) + 



^o(r)(24) 



rr/i 

j dt' In(cosht') + 5o(r) , t . (25) 



Vm2 + e2 12 J \ ^ ' ' V m? + e^E'^e 

where Sq{t) is a certain real entire function of r. There are two important limits of (j24jl . 
At t — s> oo, we can expand the integrand of (j24|) to the terms ~ 1/t. Then the imaginary 
part of the action Sj reduces to ImS'/(r) = iT'm/2eE which coincides with the corresponding 
result in the constant field. 
In the opposite limit t ^ 

Sj(r) = ^^^^ 
m 

The imaginary part of is calculated by making transformation to the Euclidean time 
tE and integrating in the interval — 7r/2 < < 7r/2, i. e. between the turning points of the 
potential. The result is 

lmSi(r) = ^ ^7rln2, t ^ 0. (26) 

m 

In the limit t ^ a^^ = m/eE, i. e. when the field E is quasi-constant we have 

2 

Tfl 

ImSir) =lmSK{r) +lmSi{T) = —j^, t ^ oo . (27) 

z \e\ hi 

Therefore we recover the previous result (jllj) . 

In the opposite limit of a short pulse the field E varies rapidly, so that t ^ a~^, i. e. the 
time t is much shorter than the inverse acceleration the particle would have obtained in a 
constant field. In that case IvuSk ^ ImS"/ as is seen from and (j2HI)- Thus, 



Im5(r) = vrmt, at < 1. (28) 

Therefore, the decay rate 

rv^„ ~ e-2--\ at<l (29) 

becomes independent of the value of the field E in agreement with the quantum mechanical 
uncertainty relation. The result (|29|) agrees with the at ^ 1 limit of the exact result 
obtained in |3Q]. Our quasi-classical derivation suggests that is a general result for 
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all electric field pulses of a typical duration t. We have verified by explicit calculations for 
several functional dependencies of E{t) that this is indeed the case. 

We can see that the distribution (f^n|) is again thermal, with an effective temperature 
of (27rt)^^. This result reflects the fact that the thermal character of the distribution is a 
general property of an accelerated reference frame, and is not an artifact of a particular 
assumption about the field which causes the acceleration. 

C. Pair production in a chromo-electric field 

The results presented above can be generalized to the case of a color particle movin g in 
a homogeneous chromoelectric field. To this end let us consider the Wong equations |37l | 
governing the classical motion of such a particle 

mx^ = gF^^^x^Ia, (30a) 



- gfabcX^A^r = 0, (30b) 

where a dot denotes the derivative with respect to the proper time and x^ = 1. The 
chromoelectric field is given by 

El ^ = d,A\ - dal + gfa,cAWc- (31) 

The homogeneous chromoelectric field can be described by the potential 

Al = -Ez5''\ K = f). (32) 

Substituting (jH^ into (j30bj) we find that the color isospin vector la precesses about the 
3-axis with = const. Therefore, ()30ap tells us that 

X = y = , m z = g E x^ I3 . (33) 

Mathematically, equations are equivalent to the equations of motion of an electric charge 
(7/3 in the electric field Ez. Therefore, all arguments given above are valid for the motion in 
a homogeneous electric field as well. Heavy quark production by Schwinger mechanism in 
strong color fields within the Color Glass Condensate picture has been previously discussed 
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The constant in time field can also be described with another choice of potentials, which 
is not related by the gauge transformation to (IH^ 139^ . However, that choice leads to the 
periodic motion of the particle along the z-axes |4y] and is of no interest for the physical 
process we will discuss in the next section. 

III. EVENT HORIZON AND THERMALIZATION IN HIGH ENERGY 

HADRONIC INTERACTIONS 

A. Unruh effect, Hagedorn temperature, and parton saturation 

We are now ready to address the case of hadron interactions at high energies, which is 
the main subject of this paper. Consider a high-energy hadron of mass m and momentum P 
which interacts with an external field (e.g., another hadron) and transforms into a hadronic 
final state of invariant mass M ^ m. This transformation is accompanied by a change in 
the longitudinal momentum 

2 ^2 



- m 



ql = VE^ - m2 - _ ^2 ^ _ ^ (34) 

and therefore by a deceleration; we assumed that the particle m is relativistic, with energy 
E ~ p. 

The probability for a transition to a state with an invariant mass M is given by 

P(M^m) = 27r|r(M^m)pp(M), (35) 

where T(M ^ m) is the transition amplitude, and p{M) is the density of hadronic final 
states^. According to the results of the previous section, we expect that under the influence 
of deceleration a which accompanies the change of momentum (jH^ . the probability |Tp will 
be determined by the Unruh effect and will be given by 

|T(M ^ m)|2 ~ exp(-27rM/a) (36) 

in the absence of any dynamical correlations; we assume M ^ m. 

To evaluate the density of states p{M), let us first use the dual resonance model (see e.g. 



4l|, |4^), in which 

p(M) ~ exp 76 , (37) 



p{M) should not be confused with the Rindler coordinate p used in other sections of the paper. 
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where b is the universal slope of the Regge trajectories, related to the string tension a by 
the relation a = 1/ (27r6). 

The unitarity dictates that the sum of the probabilities over all finite states M should 
be finite. Therefore, by converting the sum into integral over M one can see that the eqs 
and (jHTjl impose the following bound on the value of acceleration a: 

2^^^^l^7i^^''"- „ '''' 

The quantity on the r.h.s. of (jHHj) is known as the Hagedorn temperature 43| - the "limiting 
temperature of hadronic matter" derived traditionally from hadron thermodynamics. In our 
case it stems from the existence of a "limiting acceleration" oq: 

«o = y|rv^. (39) 

The meaning of the "limiting temperature" in hadron thermodynamics is well-known: 
above it, hadronic matter undergoes a phase transition into the deconfined phase, in which 
the quarks and gluons become the dynamical degrees of freedom. To establish the meaning 
of the limiting acceleration ()39|). let us consider a dissociation of the incident hadron into a 
large number n 1 of partons. In this case the phase space density ^7\\ can be evaluated 
by the saddle point method ("statistical approximation"), with the result (see e.g. 

p{M) ~ exp(/5M); (40) 

where f3~^ is determined by a typical parton momentum in the center-of-mass frame of 
the partonic configuration. When interpreted in partonic language, eq(|H7j) thus implies a 
constant value of mean transverse momentum ~ ~ 6^^^^. On the other hand, in the 
parton saturation picture, the mean transverse momentum has to be associated with the 
"saturation scale" Qs determined by the density of partons in the transverse plane within 
the wave function of the incident hadron (or a nucleus). This leads to the phase space 
density logp(M) ~ M/Qg. The unitarity condition and the formulae (jHSI) thus lead 

us to the acceleration a ^ Qg, which can exceed (j39|) . and to the conclusion that the final 
partonic states are described by a thermal distribution with the temperature T ~ Qs/{27t). 
The same result can be obtained by considering the acceleration © a = gE/m of a parton 
with off-shellness m = a/p^ — Qs in an external color field gE ^ Ql. Using Wong equations 
fj33|l we derive a = Qs and 

T-9l ^A^^ 
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It is interesting to note that to exceed the hmiting acceleration ()39p. and thus the hmiting 
Hagedorn temperature (jHHj) for the produced hadronic matter, one has to build up strong 
color fields, exceeding qEq ~ 1/6. This is achieved by parton saturation in the Color Glass 
Condensate, when gE ^ > qEq at sufficiently high energies and/or large mass numbers 
of the colliding nuclei. Parton saturation in the initial wave functions thus seems to be a 
necessary pre-requisite for the emergence of thermal deconfined partonic matter in the final 
state. We will discuss this conjecture in more detail below. 

According to arguments in the previous sections, the thermal distribution is built over 
the time period of 

ttherm ~ = ^. (42) 

As discussed above, this apparent thermalization originates from the presence of the event 
horizon in an accelerating frame: the incident hadron decelerates in an external color field, 
which causes the emergence of the causal horizon. Quantum tunneling through this event 
horizon then produces a thermal final state of partons, in complete analogy with the thermal 
character of quantum radiation from black holes. 



B. The space time picture of relativistic heavy ion collisions 

The conventional space-time picture of a relativistic heavy ion collision is depicted in the 
left panel of FigI21 The colliding heavy ions approach the interaction region along the light 
cone from x = t = —oo and x = —t = oo. The partons inside the nuclei in the spirit of the 
collinear factorization approach are assumed to have a vanishing transverse momentum kT, 
have a zero virtuality k"^ = —k^ = 0, and thus are also localized on the light cone at zizx — t. 
The collision aX x = t = produces the final state particles with transverse momenta pt 
which according to the uncertainty principle approach their mass shell at a proper time 



For further discussion it is convenient to introduce the Rindler coordinates 

t + X 



P 



X 



V 



2 



t 



X 



(43) 



The surface of a fixed proper time which is a hyperbola in Minkowski space thus represents 



a line at 



X 



= Po < in Rindler space. The Rindler coordinate 7] in high energy 
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FIG. 2: The space-time picture of a relativistic heavy ion collision. Left: heavy ions approach the 
interaction region around x = t = along the light cone from x = t = —oo and x = —t = oo. The 
collision at x = t = produces the final state particles which approach the mass shell at some proper 
time r = (t^— j;^)^/^ (or, equivalently, along the surface of Rindler space = Pq = x'^—t^ < 0). The 
produced particles are distributed in rapidity, or in the Rindler coordinate r/ = ^ In 
also shown in the left (L) and right (R) sectors are the trajectories of space-like ~ —Q^ < 
partons confined in the initial nuclear wave functions characterized by the saturation scale Qs- 
Approaching the interaction region around x = t = the partons from the colliding nuclei begin 
to interact, which leads to their deceleration with \a\ = Qs', the trajectories in the left and right 
sectors are the hyperbolae with p'^ = —p^ > 0, \a\ = Qs = p~^- For partons in the right (R) sector, 
the surface p'^ = 0, r] = +oo is the event horizon of the future; the information about the left (L) 
and the future (F) sectors is hidden from them. Quantum tunneling of partons from left (L) and 
right (R) sectors through this event horizon into the future (F) indicated by dashed arrows creates 
a thermal state of parton matter with the temperature T = Qs/2tt. 

physics is often called a space-time rapidity. 

Consider now the case when partons in the wave functions of the colliding nuclei have 
non-vanishing transverse momenta, as in the Color Glass Condensate picture where their 
transverse momenta are on the order of the saturation scale Qg. In this case the partons are 
space-like k'^ = —k^ and are located off the light cone. As the colliding nuclei approach each 
other, these partons begin to interact; note that since they are space-like, their interactions 



. Right: 
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are acausal, and are responsible for the breakdown of factorization for the parton modes 
with krp < Qs- The interactions of partons with the color fields gE ~ of another nucleus 
decelerate them, with a typical acceleration \a\ ^ Qs- The space-time trajectories of the 



partons, according to eq. (j8b|l are thus given by hyperbolae in Minkowski space, or by the 
lines with a fixed value of = — Pq = > in Rindler space. These trajectories are 
shown on the right panel of Fig|21 For partons moving in the left (L) and right (R) sectors 
of space-time with an acceleration \a\ = the light cone surfaces p"^ = 0, 1] = ±oo 
represent the event horizon of the future (F). The information from the future is hidden 
from them, and the sector F is classically disconnected from L and R. However, as discussed 
above in Section InJ the future sector F can be reached from the left L and right R sectors 
by quantum tunneling. When evaluated in the quasi-classical approximation, for a parton 
of mass m and transverse momentum the tunneling probability is 



WM ~ exp ( ^^g^-^ 1 • (44) 

Eq. ()44j) describes a thermal distribution in the "transverse mass" vtit = \/ m"^ + Pr with the 
temperature T = Qs/2-n. 

The partons produced in the sector of future F reach the mass shell at a typical proper 
time To = (— Po)^^^ = 1/Qs- Note that the Rindler trajectory p^ = — Pq of the produced 
partons in sector F is related to the trajectories of colliding partons p^ = Po the sectors 
L and R by crossing symmetry. The thermalization of the produced partons is reached over 
the time Ttherm = = 27t/Qs. 



IV. IMPLICATIONS OF THE "BLACK HOLE THERMALIZATION" FOR 
RELATIVISTIC HEAVY ION COLLISIONS 



A. Rapid thermalization 



There is an ample amount of evidence 



i],y,y,iio 



that heavy ion collisions at RHIC 



energies lead to a thermalized state of matter (see 111 and references therein). Hydrody- 

nnn 

namical models |12I. I13LI14| appear successful in describing the collective flow of the produced 
hadrons, but only if the thermalization time Ttherm (at which the hydrodynamical evolution 
begins) is short, Ttherm < 1 fm. None of the existing approaches had so far succeeded 
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in explaining such a fast thermalization, despite some promising work ranging from the 



"bottom- up" scenario 



3ased on parton saturation 

HQ, 00,0,0 



plasma instabilities 



42l 



le p 

m 



45, 4&I to the consideration of collective 



The mechanism of "black hole thermalization" proposed in this paper produces soft ther- 
mal partons with the initial temperature of T = Qs/2tt] the thermalization time in this case 
is Tfhprm = = 27r/Qo. The analysis of experimental data at RHIC in the framework of 
KLN sa..a.„n .nodel bl yield, fo. Au-Au coHi.o. the va.„e. of Q. . 0.9 . 1.5 GeV 
depending on centrality. This corresponds to the temperatures of T ^ 140 -i- 240 MeV, 
and to thermalization times of about Ttherm — 1 fm. Note that the lower bound on the 
initial temperature according to our discussion in section IIII Al is given by the Hagedorn 
temperature. Indeed, even if the saturation in the initial wave function is not reached, the 
deceleration in hadronic collisions is determined by the string dynamics; the acceleration in 
this case is determined by the string tension, related through the density of hadron states 
to the Hagedorn temperature, see 



B. Phase transitions 



The arguments presented in this paper are admittedly somewhat schematic. Nevertheless, 
we hope that our approach may provide a useful theoretical tool for understanding the 
dynamics of hadronic interactions. As an explicit example of using the geometry of Rindler 
space we will now address an important issue of phase transitions induced by relativistic 
heavy ion collisions. 

As seen from © the stronger the external field, the larger acceleration it causes since 
a oc E. The strength of the saturated color field of a nucleus increases with the energy of 
the collision squared s and the atomic number A as E (x Ql ^ A^^^s^^'^ where A ~ 0.3. 
Thus the acceleration of particles in hadron (nucleus) - hadron (nucleus) collisions grows at 
higher energies and for higher atomic masses as a power. On the other hand, Unruh effect 
relates the acceleration to the temperature of the thermal bath of particles surrounding the 
accelerated one, see (0). Therefore, by increasing the energy of the collision one increases 
the temperature of the produced thermal bath, (j4T|l . This allows tuning an order parameter 
characterizing the thermal bath of the produced particles. In particular, by varying the 
saturation scale in a system with spontaneously broken symmetry one can observe the phase 
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transition to the symmetric phase. Here we would hke to discuss an example of such a phase 
transition; to be concrete, we will concentrate on the chiral symmetry restoration. 

Rindler space is often considered in general relativity as an approximation to the 
Schwarzschild metric of a large black hole; consequently, a number of important calcula- 
tions had already been performed in this space. In this section, we will rely on these results; 
our treatment will closely follow a recent paper by Ohsaku . 

Consider the A^-fiavor Nambu-Jona-Lasinio (NJL) model: 

where is the Dirac field and A is the coupling. The gamma matrices satisfy the following 
anti-commutation relations: 

{lf,{x),-f,y{x)} = 2g^^{x). (46) 

A uniformly accelerated observer will measure the thermal spectrum of V^'s with temperature 
T = 2t{ / a. To find an effective Lagrangian of an accelerated observer we transform the NJL 
theory to the Rindler coordinates using (021) • The interval in Rindler space takes form 

ds"^ = drf - dp^ - dxl . (47) 

With the help of ()46|) one can express gamma matrices in the Rindler space through those 
in the Minkowski space. 

Introducing as usual the scalar and pseudo-scalar fields as 

(t{x) = -^ip{x)'ilj{x) , tt{x) = -^^(x)i75^(x), (48) 

and performing large- iV expansion one can calculate partition function of the problem from 
which the effective action can be read off: 



'eff 



y d^x ^^^^^ - Undet(i7'^V, -(J -i75 7r). (49) 



This action can be used to derive the gap equation 
d^k^_ p , sinh(7T 



2i\cr f d'^k_i r°° sinh(7ra;/a) , ^ i ^^2 \ /r-n\ 

a = — / / du -z^^{ {Ki^+i{a/a)y - {K,^_i{a/a)y \, (50) 
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Veff 




FIG. 3: The effective potential ()52|) of NJL model (scaled by 32it) in an accelerated frame near 
the critical point ()51() . We chose the critical value of the saturation scale Qs,c = 1 GeV which is 
equivalent to the critical temperature of Tc ~ 160 MeV. 

where = /c^ + cr^. The critical value of the acceleration Oc at which the phase transition 
occurs is determined by the gap equation (jSUj) at a = 0. One gets 

where A is the cutoff. To study the system near the phase transition a ~ Oc one expands 
(I5()|) in the vicinity of the critical point 1)511) . The effective potential reads: 

VMQs, ^) = -w^ {Qlc -Ql)a' + -^a' + ..., (52) 

where we used ()4H) . Eq. ()52j) demonstrates that chiral symmetry restoration occurs whenever 
the center of mass energy ^/s or nucleus atomic number A is high enough to ensure that the 
saturation scale Qs is larger than the critical value Qs,c = clc, see Fig. 01 

We believe that a similar treatment can be extended to other phase transitions in heavy 
ion collisions, including the deconfinement and the Ua{^) symmetry restoration. 
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C. Hadron abundances, spectra, and the HBT radii 

It has been known for a long time since the pioneering work by Fermi and Landau 
that statistical models successfully describe many features of hadronic reactions 
at high energies. Indeed, the abundances of produced hadron s by and large follow the 
statistical ones with an effective temperature T ~ 150 MeV (see js^ 58, and references 



therein); moreover, the hadron spectra at small "transverse mass" tjit = \/ vrP- + look 
approximately thermal, ~ exp(— m^'/T). 

On the other hand, it has been found that the variation of the transverse momentum 
spectra at RHIC with centrality (the "transverse flow") exhibits an approximate scaling 
pd l6l| in niT/Qs- The picture proposed in this paper provides a natural explanation to 
this fact: a thermal spectrum ()44p with T = (5s/27r obviously scales in tjit/Qsi even though 
the shape of the thermal distribution of course can be affected by the hydrodynamical 
evolution. 

Moreover, it has been found both in 125^ and in 0| that the saturation scale has to 
be "frozen" at low energies and/or peripheral AA collisions at a lower bound of about 
QgQ ^ 0.9 GeV to describe the data. In view of our present discussion, this cutoff is related 
to the critical acceleration, or the Hagedorn temperature: Qso = = '^T^Tuagedorn- This 
may explain why the elementary e^e^ and p{p)p collisions exhibit many statistical features 



63|: the "subcritical" color fields with the strength described in the dual resonance 



model by the string tension induce the acceleration a ~ Oc = ^nTHagedom- 

A well-known puzzle stemming from RHIC results is the short duration of hadron emis- 
sion as implied by the ratio of HBT radii Rside/Rout — 1- This is difficult to achieve if 
the thermalization is reached by the consecutive interactions of the produced particles. We 
speculate that our scenario with a short production time t = 2'k/Qs oi the apparently ther- 
malized particles may help to solve this problem, but our claim of course has to be verified 
by an explicit calculation. 



V. SUMMARY AND DISCUSSION 



In this paper we have proposed a new "black hole thermalization" scenario for hadron 
and heavy ion collisions. The key idea of our approach is the existence of an event horizon 
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caused by a rapid deceleration of the colliding nuclei. The apparent thermalization in this 
case is a general phenomenon known as the Hawking-Unruh effect: an observer moving with 
an acceleration a experiences the influence of a thermal bath with an effective temperature 
T = a/2TT, similar to the one present in the vicinity of a black hole horizon. In hadron 
collisions, the deceleration Oc is given by Eq. (jSH)), and appears related to the Hagedorn 
temperature: Oc = ^nTHagedorn- To exceed the Hagedorn temperature, and to induce the 
phase transition(s) to the deconfined and chirally symmetric phase, one needs a larger de- 
celeration which can be achieved if the color fields within the colliding hadrons or nuclei 
grow stronger. A mechanism leading to strong color fields in the initial wave functions is 
at present well-known: it is parton saturation in the Color Glass Condensate. In this pic- 
ture, the strength of the color-electric field is ~ Q^/g {Qs is the saturation scale, and 
g is the strong coupling), the typical acceleration is a ~ Qs, and the resulting heat bath 
temperature is T = Qs/27r ~ 200 MeV. Such a large deceleration should induce a rapid 
thermalization in nuclear collisions over the time period of r ~ 27r/Qs ~ 1 fm accompanied 
by phase transitions. 

We have considered an explicit example of chiral symmetry restoration induced by a rapid 
deceleration of the colliding nuclei. On these grounds, we have argued that parton saturation 
in the initial nuclear wave functions is a necessary pre-condition for the formation of quark- 
gluon plasma. We have discussed the possible implications of our "black hole thermalization" 
scenario for various observables in relativistic heavy ion collisions, including hadron spectra 
and abundances, and the HBT radii. 

Our discussion of the Schwinger mechanism in section |n] was aimed at the understanding 
of particle production in strong color fields in the framework of the geometrical Hawking- 
Unruh approach. We have argued that the pair production in this case can be understood as 
a quantum tunneling through the event horizon, similar to the quantum evaporation of black 
holes. In our case, the event horizon appears due to the acceleration of particles in external 
fields. Schwinger mechanism has been extensively discussed as a model for multi-particle 
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6q . It was also argued that this 



production in hadronic and nuclear interactions 
pair production can drive the system towards equilibrium through successive interactions. 
Our proposal here is quite different: we show that the spectrum of the produced particles 
is inherently thermal, and that this apparent thermalization is achieved over a very short 
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time of tform ^27l/Qs. 

Numerical simulations of the classical Yang-Mills dynamics in nuclear collisions [gtI . IgsI ] 
were found to lead to approximately thermal transverse momentum spectra of the produced 
gluons, with an effective temperature T ~ Qs/2, substantially higher than ours. In our case, 
the thermal spectrum arises due to a purely quantum effect; since the quantum behavior 
cannot be reproduced by a classical simulation, our results have to be of different physical 
origin. 

A shortcoming of the quasi-classical approach used in this paper is the necessity to con- 
sider the trajectories of particles propagating in external classical color fields. On the other 
hand, the classical color fields themselves are formed by colored particles - gluons, with 
the occupation numbers ~ l/f?^- The necessity to treat differently the "particles" and the 
"classical fields" of course is common for many kinetic approaches. In our case, it may be 
possible to overcome this difficulty by formulating the field theory in Rindler space. We 
have provided an explicit example of such an approach in section IIVBI 
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